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The  state-of-the-art  hypergol  combination  currently  used  in  the  US  for  many  space  propulsion 
applications  consists  of  monomethyl  hydrazine,  as  the  fuel,  and  nitrogen  tetroxide,  as  the  oxidizer. 
The  Air  Force  Research  Laboratory  is  developing  new  hypergolic  fuels  which  will  provide  enhanced 
performance  capabilities  as  well  as  improved  affordability  and  efficiency.  Furthermore,  handling  of 
these  new  hypergolic  fuels  is  also  expected  to  have  a  much  smaller  logistical  footprint  due  to  the  fact 
that  they  are  being  designed  to  be  environmentally  benign.  Flowever,  practical  realization  of  these 
hypergols  in  spacecraft  propulsion  systems  will  only  come  after  attaining  a  satisfactory  understanding 
of  how  to  optimize  their  combustion  characteristics  in  relevant  operating  environments.  Flere  we 
report  theoretical  results  obtained  on  the  prototypical  radical-radical  reaction:  NCb  +  N2H3,  and  the 
progress  made  towards  building  an  apparatus  consisting  of  laser  photolysis/fast  flow-tube  reactor 
coupled  to  a  mass  spectrometer  for  investigating  the  kinetics  of  this  elementary  reaction. 
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1.  Introduction 

The  hydrazine  (N2H4)  and  dinitrogen  tertroxide  (N2O4)  reaction  system  is  widely  used  as  liquid 
propellants  for  satellite  and  space  propulsion  applications.1  As  an  oxidizer,  the  N2O4  molecule  is 
known  to  have  several  inter-conversion  structural  confonners:  NO2  monomer,  symmetric  N2O4 
dimer  with  D2h  symmetry,  asymmetric  ONONO2  with  cis-  and  trans-structures,  and  the  van  der 
Waals  complex  with  CLv  symmetry."  Recently,  Lai  et  al.  theoretically  studied  the  reaction  of  N2H4 
with  the  N2O4  isomers  and  reported  that  the  reaction  of  N2H4  with  sym-N204  (D21O,  cis-ONONC>2 
(Cs),  and  NO2  has  to  overcome  an  activation  barrier  of  14.2,  10.6,  and  7.6  kcal/mol,  respectively. 
However,  N2H4  can  spontaneously  react  with  trans-0N0N02  (Cs)  and  cis-ONONCL  (Ci)  to  produce 
HONO2  +  H2NN(H)NO  with  a  rate  constant  of  4  x  10'10  cm3  molecule'1  s'1  at  250  K.  The 
intermediate  H2NN(H)NO  thus  formed  can  rapidly  fragment  to  produce  the  reactive  radical  N2H3 
with  a  rate  constant  of  1  x  107  s'1  at  1000  K,  which  is  fast  enough  to  initiate  the  hypergolic  chain 
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reactions  of  NO2  and  N2O4.  Consequently,  the  kinetics  of  N2H3  with  NO2  is  important  for 
understanding  the  mechanism  and  the  ignition  process  of  the  hydrazine  and  N2O4  system.  However, 
the  reaction  is  complicated  by  multiple  radical-radical  abstraction,  addition  and  subsequent 
dissociation  pathways.  In  this  work,  the  reaction  of  N2H3  with  NO2  was  studied  by  ab  initio  multi¬ 
reference  second-order  perturbation  theory  and  quadratic  configuration  interactions  and  coupled- 
cluster  theories.  It  was  found  that  oxidation  of  the  N2H3  radical  by  NO2  occurs  with  submerged 
barriers  with  relatively  large  exothermicities.  Consequently,  the  accuracy  of  the  value  of  the 
theoretical  rate  coefficient  computed  for  this  reaction  will  greatly  depend  on  the  level  of  accuracy 
employed  to  describe  the  electron  correlation  effects  for  both  the  outer  and  multiple  inner  transition 
states.4  We  have  accurately  characterized  the  potential  energy  surface  for  the  NO2  +  N2H3  system 
and  have  established  the  most  likely  reaction  mechanism.  The  technique  of  laser  photolysis/fast 
flow-tube  reactor  coupled  to  a  mass  spectrometer  is  also  being  utilized  to  experimentally  investigate 
the  reaction  kinetics  of  this  system. 

2.  Theoretical  and  Computational  Details 

Due  to  electron  repulsion  between  the  electronegative  oxygen  and  nitrogen  atoms  in  this  reaction 
system,  the  multi-reference  character  of  the  wavefunction  is  quite  strong.  Therefore,  we  have  used 
multi-reference  second-order  perturbation  theory  (CASPT2)  with  Dunning’s  augmented  correlation 
consistent  basis  set,  aug-cc-pVDZ  to  optimize  the  geometries  of  the  stationary  points  of  the  potential 
energy  surface  and  to  calculate  the  corresponding  ro-vibrational  frequencies.  For  the  reactants 
(N2H3  and  NO2)  and  the  products  (NNH2,  cis-NH=NH,  and  trans-NH=NH),  the  larger  basis  set,  aug- 
cc-pVTZ  was  applied,  and  all  of  the  electrons  were  included  in  the  active  space.  The  geometries  of 
the  products  (cis-HONO,  trans-HONO,  and  HN02)  were  optimized  at  the  CASPT2/aug-cc-pVDZ 
level,  with  all  the  electrons  included  in  the  active  space.  For  the  transition  states,  the  active  space 
was  selected  to  describe  adequately  the  bond  breaking/making  process  of  interest  with  reasonable 
computational  efficiency.  Specifically,  for  H-abstraction  from  N2H3  by  NO2  to  form  cis-HONO  and 
the  NNFF  diradical,  the  transition  state  was  optimized  using  an  active  space  (6e,8o)  consisting  of  a 
bonding  and  anti-bonding  pairs  of  the  imidogen-H  being  abstracted,  n  bonding  and  anti-bonding 
pairs  of  the  N-N  bond,  s  orbital  of  the  electron-deficient  N  atom  on  NNH2,  and  a  bonding  of  the  O- 
H  bond  being  formed.  For  H-abstraction  from  N2H3  by  NO2  to  form  cis-HONO  and  trans-NH=NH, 
the  transition  state  was  optimized  using  an  active  space  (6e,6o)  consisting  of  a  and  n  bonding  and 
anti-bonding  pairs  of  N-0  bond,  p  orbital  of  the  O  atom  that  is  interacting  with  the  H  being 
abstracted  and  its  corresponding  virtual  orbital.  Furthermore,  the  (6e,8o)  active  space  was  used  for 
optimization  of  the  addition  adducts  and  H-bonded  complexes. 

Higher-level  stationary  point  energies  were  obtained  from  quadratic  configuration  interactions  and 
coupled-cluster  theories  with  single  and  double  excitations,  and  correction  for  triple  excitations. 
Specifically,  the  QCISD(T)  and  CCSD(T)  calculations  employed  the  correlation-consistent, 
polarized-valence,  triplet-^  (cc-pVTZ)  and  quadruple-^  (cc-pVQZ)  basis  sets  of  Dunning.5,6  The 
energies  were  then  extrapolated  to  the  complete  basis  set  (cc-pVooZ)  limit7  by  the  asymptotic 
form.8,9  Electronic  structure  calculations  were  performed  using  the  MOLPRO10  quantum  chemistry 
package. 
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3.  Experimental  Details 


A  halo-carbon  wax  coated  Pyrex  flow-tube  reactor  which  was  100  cm  long  with  internal  diameter  of 
2.54  cm  was  used.  N2H4  vapor  and  gaseous  NO2  were  carried  into  the  flow-tube  using  excess 
nitrogen  gas.  The  N2H4  was  photolyzed  under  optically  thin  conditions  by  sending  a  193-nm  laser 
beam  along  the  axis  of  the  flow-tube  which  had  been  fitted  with  appropriate  quartz  optical  ports  at 
both  ends  to  produce  N2H3  radicals  under  excess  NO2  conditions.  The  pressure  of  the  flow-tube  was 
typically  held  at  ~  2  Torr  and  the  linear  flow  velocity  of  the  gas  mixture  adjusted  to  attain  plug-flow 
conditions  that  gave  reaction  lengths  ranging  from  40  to  70  msec.  The  reaction  kinetics  of  the  N2H3 
radicals  was  followed  in  real-time  by  sampling  the  gas  mixture  through  a  small  pin  hole  at  the  mid¬ 
point  down  the  flow-tube  axis  using  a  skimmer  that  was  attached  to  a  differentially  pumped  chamber 
containing  an  electron  impact  ionization  quadrupole  mass  spectrometer. 


4.  Results  and  Discussion 

To  our  knowledge,  experimental  and/or  theoretical  studies  on  the  kinetics  of  the  reaction  NO2  +  N2H3 
has  not  been  reported  in  the  literature.  Recently,  Liu  et  al.11  studied  the  gas-phase  chemistry  of  the 
hypergolic  system,  monomethyl  hydrazine  (MMH)  with  NO2/N2O4  at  room  temperature,  and 
proposed  that  the  oxidation  of  MMH  by  NO2/N2O4  at  atmospheric  pressure  occurs  via  two 
mechanisms:  (1)  sequential  H-abstraction  and  HONO  formation,  and  (2)  reaction  of  MMH  with 
asymmetric  ONONO2,  leading  to  the  formation  of  methyl  nitrate.  Furthermore,  they  reported  that  all 
of  the  intermediates  observed  experimentally  could  be  successfully  explained  by  the  above 
mechanisms.  Furthermore,  Catoire  et  al.  suggested  that  the  addition  products  CH3N(NC>2)NH2  and 
CH3N(ONO)NH2  are  the  major  stable  pre-ignition  products  which  account  for  the  formation  of  the 
white -yellow  condensate.13  However,  in  the  study  of  Liu  et  al.,11  these  species  were  not  identified  as 
the  major  early  species  and  the  only  major  product  observed  is  methyldiazene,  CH3N=NH.  This 
indicates  that  the  direct  H-abstraction  reaction  from  the  methylhydrazyl  radical  by  NO2  leading  to 
the  fonnation  of  CH3N=NH  might  be  a  more  dominant  path  compared  to  the  radical-radical  addition 
path  followed  by  subsequent  decomposition  of  addition  products.  Consequently,  this  would  suggest 
that  this  might  also  be  the  case  for  the  analog  reaction  between  NO2  and  N2H3.  In  this  work,  the 
reaction  pathways  for  both  abstraction  and  addition  in  NO2  +  N2H3  were  theoretically  explored  at  the 
CCSD(T)/cc-pVooZ//CASP2/aug-cc-pVDZ  level. 

(I).  Direct  H-abstraction  Reactions 

The  ground  state  of  NO2  (“Ai)  has  C2V  symmetry.  For  direct  abstraction  of  hydrogen  from  N2H3, 
hydrogen  bonding  plays  an  important  role  in  the  transition  state  geometries.  It  was  determined  that 
direct  abstraction  occurs  via  four  different  transition  states,  all  with  negative  barriers  (see  Figure  1), 
to  fonn  the  following  products: 


NO2  +  N2H3  ->  trans-NH=NH  +  cis-HONO  (Rl) 

->  NNH2  +  cis-HONO  (R2) 

->  NNH2  +  trans-HONO  (R3) 

->  NNH2  +  HN02  (R4) 
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Unit:  kcal/mol 


CASPT2/aug-cc-pVDZ 

RCCSD(T)/CBS//CASPT2/aug-cc-pVDZ 

Figure  1:  Potential  energy  surface  of  NO2  +  N2H3  (abstraction)  calculated  at  the  RCCSD(T)/cc- 
pVooZ//CASPT2/aug-cc-pVDZ  level.  Energies  are  zero-point  corrected  and  relative  to  that  of 
entrance  channel  at  0  K. 

The  most  exothermic  abstraction  channel  is  reaction  (Rl),  which  has  a  submerged  barrier  of  -5.70 
kcal/mol  via  TSAI,  forming  an  H-bonded  product  complex  (PHBC1)  with  an  energy  of  -39.57 
kcal/mol.  It  then  dissociates  to  the  products  trans-NH=NH  and  cis-HONO  with  an  overall  reaction 
energy  of  -34.55  kcal/mol.  The  abstraction  of  imidogen-H  (R2)  was  found  to  have  the  lowest  energy 
barrier  of  -12.88  kcal/mol  via  TSA2,  which  is  due  to  electron  donation  from  the  N  atom  to  the 
forming  O — H  bond.  The  reaction  forms  an  H-bonded  product  complex  (PHBC2)  with  an  energy  of 
-22.76  kcal/mol,  which  then  dissociates  to  the  products  NNH2  plus  cis-HONO  with  an  overall 
reaction  energy  of -10.94  kcal/mol.  Reaction  (R3)  has  an  energy  barrier  of  -3.17  kcal/mol  via  TSA3 
and  it  forms  an  H-bonded  product  complex  (PHBC3)  with  energy  of  -19.39  kcal/mol,  which  then 
dissociates  to  the  products  NNH2  plus  trans-HONO  with  overall  reaction  energy  of -10.84  kcal/mol. 
The  abstraction  reaction  by  the  less  nucleophilic  nitrogen  of  NO2  was  found  to  be  almost  thermal- 
neutral,  with  energy  of  reaction  of  -2.73  kcal/mol.  Even  though  the  nitrogen  is  less  nucleophilic  than 
the  oxygen  in  NO2,  the  H-abstraction  by  the  nitrogen  of  N02  still  has  a  negative  energy  barrier  of  - 
5.37  kcal/mol  via  TSA4  due  to  electron  donation  from  the  N  atom  to  the  forming  N — H  bond.  The 
formed  H-bonded  product  complex  (PHBC4)  from  this  channel  has  an  energy  of  -14.15  kcal/mol. 
Overall,  the  four  H-bonded  product  complexes  are  5  to  12  kcal/mol  more  stable  than  the 
corresponding  dissociation  products,  and  it  is  noted  that,  at  low  and  intermediate  pressures,  cis- 
HONO,  trans-HONO,  trans-NH=NH,  and  NNH^  are  the  favorable  products. 
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(II).  N-N  Addition  Reactions 


It  was  found  that  both  the  nitrogen  and  the  oxygen  of  the  NO2  can  be  added  to  the  N2H3  radical.  For 
the  case  of  nitrogen  addition  to  the  radical  site  in  N2H3  (N-N  addition),  as  shown  in  Figure  2,  two 
adducts  are  formed,  which  are  optical  isomers  with  energy  of  -37.9  kcal/mol. 


CASPT2/aug-cc-pVDZ 

RCCSD(T)/CBS//CASPT2(8e,6o)/aug-cc-pVDZ 

RQCISD(T)/CBS//CASPT2(8e,6o)/aug-cc-pVDZ 

Figure  2:  Potential  Energy  surface  of  NO2  +  N2H3  (addition)  calculated  at  the  RCCSD(T)/cc- 
pVooZ//CASPT2/aug-cc-pVDZ  level.  Energies  are  zero-point  corrected  and  relative  to  that  of 
entrance  channel  at  0  K. 

In  Figure  2,  adduct  Addl  is  fonned  with  the  imidogen-H  in  the  front  of  the  pseudo  N-N-N-0  plane, 
and  Add2  is  formed  with  the  imidogen-H  in  the  back  of  the  pseudo  N-N-N-0  plane.  By  intrinsic 
reaction  coordinate  (IRC)  analysis,  it  was  found  that  the  two  isomeric  adducts  are  formed  via  one 
saddle  point,  TS1,  which  has  a  planar  structure  with  an  energy  barrier  of  -11.73  kcal/mol. 
Furthermore,  no  pre-reactive  H-bonded  complex  is  formed  before  the  saddle  point  due  to  the 
presence  of  the  strong  nucleophilic  O  atom  in  NO2.  It  was  also  found  that  each  of  the  optical  isomer 
adducts,  NH2NHNO2  have  distinct  transition  states  corresponding  to  intramolecular  H-abstraction 
from  the  NH2  group  to  form  trans-  and  cis-NH=NH,  and  trans-HONO  products.  The  Addl  adduct 
undergoes  intramolecular  H-abstraction  via  TS2  with  an  energy  barrier  of  -9.38  kcal/mol  to  form 
trans-NH=NH  +  trans-HONO  products,  and  it  also  undergoes  an  intramolecular  H-abstraction  via 
TS3  with  an  energy  barrier  of  -2.02  kcal/mol  to  form  cis-NH=NH  +  trans-HONO  products. 
Similarly,  the  Add2  adduct  undergoes  intramolecular  H-abstraction  via  TS4  with  an  energy  barrier  of 
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-4.51  kcal/mol  to  form  cis-NH=NH  +  trans-HONO  products,  and  it  also  undergoes  intramolecular  Id- 
abstraction  via  TS5  with  an  energy  barrier  of  -4.90  kcal/mol  to  form  trans-NH=NH  +  trans-HONO 
products.  For  the  intramolecular  abstraction  of  imidogen-H,  both  adducts  share  the  same  four- 
member-ring  transition  state  TS6  which  has  an  energy  barrier  nearly  identical  to  that  of  the  entrance 
channel.  This  reaction  channel  forms  NH2NNO2H  as  the  intermediate  with  energy  of  -33.81 
kcal/mol,  which  is  stabilized  by  two  internal  H-bonds.  The  predominant  channel  in  Figure  2  is  the 
dissociation  of  Addl  via  TS2,  with  an  energy  barrier  of  -9.38  kcal/mol,  which  forms  trans-HONO  + 
trans-NH=NH  as  the  products. 

For  the  case  of  oxygen  addition  to  the  N2H3  radical,  two  adducts  NH(ONO)NH2  and  NHNH2(ONO) 
are  formed,  which  undergo  further  decomposition  reactions.  The  potential  energy  profiles  for  oxygen 
addition  pathways  are  currently  under  investigation. 

(III).  Experimental  Reaction  Kinetics 

Preliminary  experimentations  are  currently  underway  to  measure  the  rate  of  HONO  production 
under  pseudo-first-order  conditions  of  [N2H3]  in  [NO2]. 


5.  Conclusions 

The  potential  energy  profile  of  the  NO2  +  N2H3  reaction  was  investigated  by  ab  initio  multi-reference 
second-order  perturbation  theory  and  quadratic  configuration  interactions  and  coupled-cluster 
theories.  For  direct  abstraction  reactions,  it  was  found  that  the  most  exothennic  channel  involves  the 
formation  of  the  trans-NH=NH  plus  cis-HONO  products  with  reaction  exothermicity  of  34.55 
kcal/mol.  The  abstraction  of  the  imidogen-H  of  N2H3  has  the  lowest  energy  barrier  of  -12.88 
kcal/mol,  forming  the  NNH2  plus  trans-HONO  products.  For  the  addition  of  NO2  to  the  N2H3 
radical,  it  was  found  that  the  reaction  proceeds  via  a  complex  mechanism.  Specifically,  both  the 
nitrogen  and  the  oxygen  atoms  of  the  NO2  can  add  to  the  N2H3  radical.  The  N-N  addition  adducts 
were  found  to  have  an  energy  of  -37.7  kcal/mol,  which  undergo  further  dissociation  reactions  by 
intramolecular  H  transfer.  The  predominant  channel  for  the  dissociation  of  the  N-N  addition  adduct 
has  an  energy  barrier  of  -9.5  kcal/mol,  forming  the  trans-HONO  plus  trans-NH=NH  products. 

The  major  products  of  the  reaction  system  are  the  trans-  and  cis-  forms  of  HONO  and  NH=NH,  and 
NNH2.  The  negative  energy  barriers  and  large  exothennicities  of  the  reaction  channels  implicate 
significant  importance  in  the  early  stages  of  hypergolic  ignition  of  hydrazine  and  its  derivatives. 
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